Zero-energy states in triangular and trapezoidal graphene structures 
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We derive analytical solutions for the zero-energy states of degenerate shell obtained as a singu- 
lar eigenevalue problem found in tight-binding (TB) Hamiltonian of triangular graphene quantum 
dots with zigzag edges. These analytical solutions are in agreement with previous TB and density 
functional theory (DFT) results for small graphene triangles and extend to arbitrary size. We also 
generalize these solutions to trapezoidal structure which allow us to study bowtie graphene devices. 



Low dimensional graphene nanostructures are promis- 
ing candidates as building blocks for future nanoelec- 
tronic applications due to their band gaps and magnetic 
properties tunable with size and shape []], 0, 0, 0,H, @ • 
Remarkable progress has been made in cutting graphene 
sheets into nanostructures with desired shape and size, 
significantly influencing their properties [jj @, 0|. In 
particular, the existence of a band of degenerate states 
near Fermi level localized at the edges in zigzag ribbons 
@, i, [Tol and triangular dots [U, El El EI El EE EE] 
was predicted by tight-binding model and confirmed by 
density functional theory calculations. These zero-energy 
edge states play important role due to their large contri- 
bution to the density of states [HI, E3, EH • In triangular 
graphene quantum dots, numerical results show that the 
degeneracy of the band of zero-energy states is propor- 
tional to the edge size and can be made macroscopic. 
This opens up the possibility to design a strongly corre- 
lated electronic system as a function of filling of the shell, 
in analogy to the fractional quantum Hall effect [17]. 

While the existence of zero-energy states was pre- 
dicted analytically for the zigzag ribbons [8], for trian- 
gular structures, the analysis of the zero-energy states 
was limited to numerical techniques such as tight-binding 
and density functional theory for specific and small sizes 
of quantum dots. A size-independent general analyt- 
ical analysis is therefore desirable. In this work, we 
present analytical solutions to zero-energy edge states in 
graphene triangles with zigzag edges. We also show how 
the results can be generalized to the trapezoidal struc- 
tures and applied to the bowtie structures [l9[. Our 
method allows the prediction of the number of zero- 
energy states as a function of the size in all triangular, 
trapezoidal and bowtie structures. 

Our starting point is the nearest-neighbour tight bind- 
ing model. It has been successfully used to describe 
graphene lattice [20[ and applied to other graphene mate- 
rials such as nanotubes, nanoribbons and quantum dots 
@, i, E3, EH, El EI El, EH . The Hamiltonian is written 



as 



where t is hopping integral, a\ and ai are creation and 
annihilation operators on a site i respectively, and (i, j) 
indicate summation over nearest-neighbours. It is impor- 
tant to distinguish between two types of atoms which ap- 
pear in the unit cell of the honeycomb lattice of graphene 
sheet. For triangular structures, these atoms form two 
non-equivalent sublattices (A and B) and they are indi- 
cated by red and blue circles of the graphene triangle in 
Fig. [TJ Our goal is to find zero-energy solutions to the 
singular eigenvalue problem, 

= 0. 

In this case there is no coupling between two sublattices 
and the solutions can be written separately for A-type 
and 5-type atoms as = with \i = A,B. The 

coefficients q obey 



0, 



(1) 



(id) 



where the summation is over i-th nearest neighbours 
of an atom j. In other words, the sum of coefficients 
around each site must vanish [8]. Let's first focus on 
the sublattice labelled by A, represented by red atoms 
in Fig. [TJ We label each atom by two integer num- 
bers n and m (with < n, m < TV + 1, where N is the 
number of A-type atoms on the one edge). The dash 
lines and open circles indicate auxiliary atoms which will 
later help to introduce boundary conditions. We will 
now show that coefficients c n ^ m for all atoms in the tri- 
angle can be expressed as a linear combination of coef- 
ficients corresponding to atoms on one edge, i.e. c n ^. 
Starting from the first row and using Eq. (pQ), we can 
obtain all coefficients corresponding to atoms in the sec- 
ond row. For the first two coefficients from the left we 
obtain c ,i = -(c ,o + ci )0 ) and c M = -(ci j0 + c 2 ,o). 
These coefficients are just equal to the sum of two up- 
per lying coefficients with the minus sign. In analogy, 
we can write expressions for all coefficients in the sec- 
ond row. In the next step, coefficients in the third row 
are expressed as a sum of two coefficients in the second 
row. For first coefficient from the left in the third row we 
obtain c ,2 = -(co,o + = (c ,o + 2ci >0 + c 2 ,o). The 
second and third ones will have similar form. By going 



FIG. 1: (Color online) Triangular zigzag graphene structure 
with N = 3 atoms on the one edge. Under each A-type 
atom (indicated by red colour) are corresponding coefficients. 
Dash lines and open circles indicate auxiliary A-type atoms 
in the three corners, which will help to introduce three bound- 
ary conditions. For zero-energy states all coefficients can be 
expressed as superpositions of coefficients corresponding to 
atoms from the one edge (upper row of atoms in our case). 



FIG. 2: (Color online) Triangular zigzag graphene structure 
from Fig. 1. Under each B-type atom (indicated by blue 
colour) are corresponding coefficients. For convenience, we 
only left coefficients corresponding to auxiliary B-type atoms 
on the right edge. For zero energy states coefficient from up- 
per left corner (co,o) determine all other coefficients in the 
structure. Introducing three boundary conditions from auxil- 
iary atoms we obtain only trivial solution; zero-energy states 
consist of only A-type atoms. 



down rows one by one, we can obtain all coefficients in 
the structure regardless of the size of the triangle. Similar 
to the construction of Pascale triangle [22], these coeffi- 
cients can be written in a suitable form using binomial 
coefficients 



(-ir£(I) c «+M- 



(2) 



Here, it is important to emphasize that the only unknown 
are the N + 2 coefficients (c n? o's) from the first row; the 
rest are expressed as their superpositions, as it is seen 
from Eq. ([2]). In addition, we must use the boundary 
conditions: the construction of the triangle requires van- 
ishing of the coefficients corresponding to auxiliary atoms 
in each corner (Fig. 1). This gives three boundary con- 
ditions (co,o = c n,o = Co, at = 0), reducing the number of 
independent coefficients to TV — 1. 

The same analysis can be done for B-type atoms in- 
dicated by blue colour. In this case, it is convenient to 
include some of boundary conditions at the beginning 
as shown in Fig. 2, where we only keep coefficients be- 
longing to auxiliary atoms on the right edge. As a con- 
sequence, the coefficient co,o determines all other coef- 
ficients in the triangle. Since there are three auxiliary 
atoms (equivalently three boundary conditions) but only 
one independent coefficient, we can not obtain any non- 
trivial solution. Hence, zero-energy states can only con- 
sist of coefficients of one type atoms - these lying on the 



edges. Now we can write general form for the eigenvec- 
tors for zero-energy states in the triangle 
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where TV is the number of atoms on the one edge and m 
is p z orbital on A-type site (n, m). In this expression the 
only N — 1 coefficients corresponding to atoms from the 
first row are independent. Thus, we can construct N — 1 
linear independent eigenvectors which span the subspace 
with zero-energy states. This is in agreement with Ref. 
P3] ~ the number of zero-energy states in the triangle is 
N — 1, where TV is the number of atoms on one edge. 

Using the Eq. §3§ we can then construct an orthonor- 
mal basis for zero-energy states. First, we make a choice 
for the N — 1 independent coefficients c n? o, from which 
we obtain N — 1 linear independent vectors, for instance 
by choosing only one nonzero coefficient for all N — 1 
collections, different one for each eigenvector. Resulting 
eigenvectors can then be orthogonalized using standard 
Gram-Schmidt process. The last step is the normalisa- 
tion K norm of the eigenvectors, using expression 
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The method for obtaining zero-energy eigenfunction 
coefficients for the triangular structures can also be ap- 
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FIG. 3: (Color online) Single particle spectrum from tight- 
binding calculations for a) 46 atoms triangle and b) 38 atoms 
trapezoid (with N row — 4) . There are four zero-energy states 
in the triangle (number of atoms on the one edge N = 5) and 
two zero-energy states in the trapezoid in agreement with our 
analysis. Changing number of rows in the trapezoid we can 
control number of zero-energy states. 



plied to trapezoidal structures (inset of Fig. 3(b)). As 
explained above, the value of the coefficients for atoms 
in a given row is sufficient to determine the coefficients 
for atoms in the lower lying row. If we stop this process 
of going down the ladder one by one at any row, we then 
obtain a trapezoidal structure. The Eq. (j3]) takes the 
following form 
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where M = min(N + 1 — 71, N row — 1) and N row is the 
number of rows in the structure (see Fig. 09(b)). In this 
case the last row contains N — N row + 2 auxiliary atoms 
which increases the number of boundary conditions. The 









FIG. 4: (Color online) Two trapezoids (left) and bowtie 
structure (right). Each of two trapezoids has one zero-energy 
state, consisting of only A-type atoms (indicated by red 
colour) for the upper trapezoid, and consisting of only B— 
type atoms (indicated by blue colour) for the lower trapezoid. 
Connecting these two systems does not affect the zero-energy 
solutions since coefficients belonging to connecting atoms are 
zeros (the four nearest atoms to the dash line). The bowtie 
structure on the right has two zero-energy states: one which 
completely lies in upper part and consists of A-type atoms, 
and second one lies in lower part and consists of B-type 
atoms. 



number of zero-energy states is then given by N row — 2 
(for N row > 1). Here we note that similar to the triangle, 
zero-energy states consist of only one type of atoms; the 
only difference is increased number of boundary condi- 
tions. In Fig. 09(a) we show tight-binding single particle 
states for triangle with N = 5 atoms on one edge. As 
expected, there are four zero-energy states. For compari- 
son, in Fig. 09(b) we show single particle states for trape- 
zoid with the same number of atoms in a first row. Here, 
there are only two zero-energy states in agreement with 
our analysis - increasing number of boundary conditions 
decrease number of zero-energy states. We note that the 
structure which consists of only two rows (the single chain 
of benzene rings, called acene) does not have zero-energy 
states while the triangular structure with N atoms on the 
one edge has maximal number of zero-energy states equal 
to N — 1. All intermediate structures (trapezoidal struc- 
tures) have number of zero-energy states in the range 
between 1 and N — 2, depending on the number of rows. 

Finally we note that the solutions of Eq. (|4]) can also 
be applied to bowtie structures [19[ . These can be treated 
as two trapezoidal structures connected by their shorter 
base, shown on Fig. [H It is important to emphasize that 
the upper trapezoid has one zero-energy state which con- 
sists of A-type atoms (red atoms) while lower trapezoid 
has one zero-energy state which consists of B-type atoms 
(blue atoms). Connecting these two systems does not af- 



4 



feet the zero-energy solutions since coefficients belonging 
to connecting atoms are zeros. Using zero-energy eigen- 
vectors for trapezoids, Eq. we obtain expressions for 
two groups of zero-energy states in the bowtie structures 
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for upper trapezoid, where A indicate A-type atoms from 
upper part and 
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for lower one, where B' indicate 5-type atoms from lower 
part. Two parts of the bowtie structure are separated by 
the dash line in Fig. [H Coefficients c n ,o( c n o) correspond 
to N A-type (TV'.B-type) atoms from the highest (lowest) 
row in the bowtie structure from Fig. [U Note that it is 
possible to use Eqs. ([5j) and (|6]) to asymmetric bowtie 
structures consisting of two different trapezoids (N ^ 
N r ). 

In summary, we derived here analytical expression for 
zero-energy states in triangular and trapezoidal graphene 
quantum dot structures. Our method allows prediction 
of the number of zero-energy states in quantum dots of 
arbitrary size which can be understood in terms of a 
competition between the number of independent coeffi- 
cients and the number of auxiliary atoms (the number of 
boundary conditions). We also showed that the number 
of zero-energy states can be controlled by changing the 
number of rows in the trapezoidal structures but does not 
depend on the number of atoms in the base of the trape- 
zoid. Finally, we applied our results to bowtie structures 
and showed that two independent groups of zero-energy 
states coexist in these systems. 
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